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THEORETICAL STUDIES OF SOLAR LASERS AND CONVERTERS 

ABSTRACT. 

The previously constructed one dimensional model for the simulated operation of an 
iodine laser assumed that the perflouroalkyl iodide gas n - C 3 F 7 / was incompressible. 
The present study removes this simplifying assumption and considers n - C 3 F 7 / as a 
compressible fluid. 

Introduction. 

The iodine laser being modeled uses the gas n - C3F7I flowing in a cylindrical tube. 
As a first approximation to the chemical kinetics and flow it was assumed that the gas 
was incompressible and operating at a constant temperature. To remove these simplifying 
assumptions we consider the effects of fluid density variation as a function of both pressure 
and temperature in the flow. In particular we have added to our previous model: 

(i) An equation of state. 

We assume the equation of state of an ideal gas and write 

P = 1 iRT (1) 

where P is the gas pressure [Pa], rj is the gas density [Kg/m 3 ], T is the temperature [°/f], 
and R = R* /296 [J/Kg°K\, and R* is the gas constant, R* = 8.317 [J/mole °K]. For 
n - C3F7J we use the value of 296 Kg/mole to convert R* to R. 

(ii) Continuity equation. 

The continuity equation expressing the conservation of mass flow can be represented 

+ *«(»?) = 0 ( 2 ) 

where rj is the fluid density [Kg/m 3 ], and V is the fluid velocity [mjs\. For steady state 
conditions and a flow in the axial direction we let V - uh and reduce the continuity 
equation to the form 

= 0 (3) 

An integration of this equation gives 

= ci (4) 


where ci is a constant of integration. 
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(iii) Momentum equation. 


The momentum equation for a control volume having a mass rjdr, where dr is an 
element of volume, is given by 


M = JJj tfrjdr 

( 5 ) 

Prom Newton’s second law we have 


'-¥-&///** 

(6) 

where D/Di is the material derivative. We have 


<*)+§-, flj Vi* 

( 7 ) 

where da is an element of surface area. The surface integral term in 

(7) represents the efflux 


of momentum through the control volume and the volume integral term in (7) represents 
the change in momentum inside the control volume. Using the Gauss-divergence theorem 
we can change the surface integral to a volume integral and write 

///<* ■ da) = jjj[v r,99} dr ( 8 ) 

Then the momentum equation becomes 

= fff + ffjiwfyfr ( 9 ) 

where vV is the dyadic u*kk for 9 = uk, and 

V • ^ (fjw 3 ) k (10) 

We also make the substitution / = /// / dr where / is the average force per unit 
volume 

as this allows the momentum equation to be written in the form 

/= III fdT = T? = HI [ Tt + £<*"*>*] iT (11) 

Neglecting viscosity and using / = -VP as the average force per unit volume which 
is due to the fluid pressure P , we have 

dP d , . , d . 2 . 

= + Tz^ (12) 


From equation (4) we use rjut = ci and examine the steady state form of the equation (12) 

to obtain 

dP dw 

dz ° l dz 


(13) 

Integrating this equation gives 


P + ciw = ca 

(14) 

where is a constant of integration. 
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(hr) Energy Equation. 

In terms of the specific enthalpy h per nnit mass, the energy equation for the fluid 
flow is __ 

where P is the pressure, T is the temperature, k is the thermal conductivity and q = g(z) 
is the radiation heat flux. In one dimension this can be written 

dh , dh dP dP . <PT , _ 

"Si + ' lu di = -di +u 7I +K d^ + l < (“) 

The specific enthalpy h can be expressed in terms of the specific heat at constant 
pressure and 

h = C p T (17) 


Using the relation 

C* p -C: = R> (18) 

where R* — 8.317 [Jjmole - 0 K] is the universal gas constant. 

The specific heat at constant pressure C* can be expressed in terms of the specific 
heat at constant volume C*. We use the numerical data 


Temperature °K 

Cl for CiF 7 I\Jlmole°K\ 

298.15 

146 

400.00 

169 

500.00 

186 


(19) 


and assume that can be obtained from the relation 

= (%(T) = a* exp(ft [T - 300)), 298.15 < T < 500 (20) 

By applying a least square’s fit of the data with the exponential curve in (20) we obtain 
the best fit parameter values of 

a, = 147.23 = 0.0012 (21) 

This gives the steady state energy equation 

where C V (T) is the correct dimensional form of C*. 

Having neglected the viscosity terms, we also neglect the effect of thermal conductivity 
and integrate the equation 

Cj (C* (T) + R\ dT + uci dw - q(z)dz = 0 (23) 

to obtain 

c, I C,(T)dT + c t B(T - 300) + ^u 1 - Q(.) = c, (24) 

where 

Q{*) = I <l{ z ) dz 

and C 4 is a constant of integration. 
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Implementation of Model. 

In order to compare and verify the computer model with laboratory data it is assumed 
that we are given the boundary conditions as z — 0 and z — L. We use the subscript 0 to 
denote values at z = 0 and subscript L to denote the values at z = L. By setting z = 0 we 
obtain: 

/rom(4) Ci = »7oWo 
from (14) 02 = Po + ciWo 

from (24) we have 

c t = c, Ji(r 0 - 300) + j C,(T)dT + ^.«g 

Assuming that g = q(z) can be represented by a constant, at z = L we obtain from the 
equation (24) the relation 

c, R)T L - 300) + Jc,{T)dT +^ L -cA 

Knowing the values of the above constants the temperature, density, pressure and flow 
rate as a function of distance z can be calculated as follows: 

From the equation of state we can write 

P = rjRT 

and multiplication by u gives 

Pw = rjwRT = ciRT 

Using the relation (14) in the above equation gives 

(ca - ciw)w = c\RT 



C2U) — CiW 2 


We also use the scale factors 

molt,K! “ d SFq ° l, ~7ffl 106) W”* , )/ r ‘ 

where applicable to insure all equations balance dimensionally. 

For 

0 < Z < 2 t 

\go*i, Zi <z<L 

We may then solve for w = u(z) such that the energy equation 


F(w) = e t j C,{T) dT + CiP(r - 300) + ^-w 3 - Q(z) -c 4 = 0 (28) 

is satisfied where T is determined by the relation (25). Since this is a nonlinear equation, 
the method of interval halving is used. The value of w is changed until P(w) changes sign 
and then interval halving is used until the desired solution accurracy is obtained. The 
temperature T is then determined from (25), the density is obtained from the relation (4) 
and the pressure is obtained from the relation (1). We then have the flow rate, temperature, 
pressure and density as functions of distance z. 


4 



Reaction Rates. 

The reaction rate coefficients 


ki } hi, £ 3 , £ 43 , £ 5 , kj, kg } 

Cl»C3|C3jC4,C6>^» 

Ql>Ql,Q3,Q*,Q6 

which occur in the chemical kinetics associated with the iodine laser can be assumed 
to vary with temperature. Various assumptions can be made concerning the effects of 
temperature on these rate coefficients. Various Arrenius type temperature variation effects 
can be modeled using one or more of the following expressions: 

K=k 0 S T ~ T '>) 

K=ko^ T 

K=k 0 T~* r 

K=kokT- N eV T 

K=hoT- N e/*l Tm 

Since the coefficients in these expressions are not known, as a first approximation to 
temperature on the rate coefficients, we will keep most of the rate coefficients constant. 

Incompressible flow model. 

The incompressible flow model was presented in the previous progress report (reference 
2). Nominal results from that model are illustrated in the figures 1 through $ for various 
values of the flow rate uj. Values of w used in these graphs are: 

w = 750,, 1750., 2750., 5500. cm/s 

while other parameters were held constant. These other parameters were: 

P = 23.5 Torr, 

R\ = .95848 reflectivity at left end of the tube, 

R? = .81634 reflectivity at right end of the tube, 
concentration = 1267 (peak concentration in solar constants.) 

Both the incompressible and compressible flow models are being modified to include 
temperature effects. Also an analysis by R.C. Costen of the pumping photon density is 
being studied and incorporated into the model. This analysis is given in the Appendix A 
along with figures of the results therein. 
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IhFOPLT 2. 

Figure 2. Log (I 2 ) vs axial distance. 


7 



LOG Cl*: 



INFOPLT 3. 

Figure 3. Log (I*) vs axial distance. 
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INFOPLT 4. 

Figure 4. Log (I) vs axial distance. 
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Appendix A 


Pumping Photon Density within Absorbing Lasant Gas 
by R.G. Costen 

Assumptions: There is a 2-D cylindrical geometry which is treated as a black body 
cavity. 

A. Consider the relation of 2-D surface brightness ^(W/m 2 ) to the 2-D surface 
flux F^W/m 2 ). Consider a one meter depth of a cylindrical cavity of arbitrary shape as 
illustrated in the figure. A.l 



We have 


Fa = 



Bq cos 9 d$ = 2 $2 


W per m of perimeter entering the cavity per m of depth. 
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B. Determination of the pumping photon density at an interior point resulting from 
a length dL of perimeter. 



Figure A.2 Density at an interior point 
The power in dB due to a length dL is given by 

Bi<x&9 dLdB 

W in d$ per unit depth. This power (J/s) is due to the photons moving at the velocity of 
light c (m/s) The energy density in d0 per m of depth is 

— cos 9dLdB (J/m) 

Using dB — and noting that the photon energy density dp? due to dL is given by 

dpa = ^ cos B dL (J/m 2 ) 
per m depth. For an absorbing gas we have 

dp 2 = ^-exp(-j) cos BdL (J/m 2 ) 

per m depth and where 5, (m) is the absorption length. 


13 



C. 2-D Circular cylinder with absorbing gas and uniform incident flux F a . 



At the point P we have due to dL that 

dp = ^-e“ * cos 6 dL {Jfv i 2 ) 

lot 

per m of depth. Using 

Si Urn# 

r T 

we can write 

°r 

P=J£ j Q = e~idi> {Jim 2 ) 

per m of depth. Using the law of cosines 


r = -ft cos ^ sin 2 i> 

we have % 

P=^-j o exP ^(ricosV»- ^/rg-rfsin 3 ^ # 

Changing fj to r the pumping photon density is given by: 

p(r|r 0> 5) = Y j[ «P cos^ - yVg-r?ain 3 ^ 
in J/m a per m of depth. 
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Figure A.4 Geometry 


D. Average pumping photon density 

The average pumping photon density is given by 


<p(ro,6)>= ^ J'° dr2nrp(r\ro,S) 

Letting s = r/6 where s 0 = r 0 /8 is the laser tube radius per absorption length, we have 

p(r|r 0 , 5) = — f exp(® cos ^ - \l*o~ & sin 2 $) dtf> 
c Jo v 

= p(s, So) 
and 

c />(ro,S) > = ^ 2* (j) *•,*.) 

2 r*° 

= 3 / da ap(s, so) 
a 0 Jo 

-< p(«o) > 

We define and compute 

f/(s, 3 0 ) = ^ J exp(s cos ^ - \]%- s a sin 3 i>) d\f> 

with 

2 f >0 

<A*o )>-Zsl *V(«j«o) 

»o Jo 

which is the average pumping photon density. 

Then 


< p(s 0 ) > — da «/(«, so) 

= < /(*o) > 
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Note: 


lim /(s,«o) = 1 
0-+ 00 


lim < /(so) > = 1 

9 0 


lim p(8,3o) 

0-+CO 

lim < p(a 0 ) >= — 


ff-Fa 

c 


and at the center s = 0 we have 

p'{ 0, So) = — f e - * 0 # = e~* 0 
* Jo 


and on the wall s = «o we have 


/^(so,*)) = \ [ e° # = 1 
* yo 
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Figure A. 5 <p ' > vs Sq. 
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Figure A. 6(a) p' vs S/S 
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Fiqure A. 7 <p 1 > vs Oil semi-log scale. 
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